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Abstract 


We compute the two-loop minimal form factors of all operators in the SU(2) sector of 
planar V = 4 SYM theory via on-shell unitarity methods. From the UV divergence of 
this result, we obtain the two-loop dilatation operator in this sector. Furthermore, we 
calculate the corresponding finite remainder functions. Since the operators break the 
supersymmetry, the remainder functions do not have the property of uniform transcen- 
dentality. However, the leading transcendentality part turns out to be universal and 
is identical to the corresponding BPS expression. The remainder functions are shown 
to satisfy linear relations which can be explained by Ward identities of form factors 
following from R-symmetry. 
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1 Introduction 


In the last years, form factors in M = 4 super Yang-Mills (SYM) theory have received 
increasing attention, both at weak coupling [1-17] and at strong coupling [18-20]. Con- 
taining both on-shell states and local composite operators, form factors provide a useful 
bridge between the purely on-shell amplitudes and the off-shell world of correlation func- 
tions. In particular, powerful computational methods developed in the context of scattering 
amplitudes can be applied to form factors and to other important physical quantities via 
form factors, such as the spectrum of anomalous dimensions of composite operators and 
their correlation functions. The form factor Fo is defined as the matrix element of a given 
local operator O(x) between the vacuum |0) and an on-shell n-particle state (1,..., n], i.e. 


Fo(1,...,2;q) = fisen... nowo. (1.1) 


A special class of form factors are the so-called minimal form factors, which contain as 
many external fields n as there are fields in the operator, and which will be of particular 
interest for this paper. 

Understanding the connection between form factors and the spectral problem of N = 4 
SYM theory was recently pushed forward in [16, 17]. In [16], form factors for generic oper- 
ators were investigated. In particular, it was shown that the complete one-loop dilatation 
operator [21] can be derived using one-loop minimal form factors, which explains the re- 
lation between the one-loop dilatation operator and the four-point scattering amplitude 


derived from symmetry in [22].! In [17], it was demonstrated that form factors can also 
be used to calculate anomalous dimensions at two-loop order by investigating the Konishi 
primary operator. In these studies, on-shell amplitude techniques have played a major 
role, in particular the (generalised) unitarity method [23-25]. In order to treat general 
operators, however, an extension of this method is required [17]. 

Interesting on-shell approaches towards the computation of correlation functions and 
the dilatation operator were also applied in the following works: see [10, 26] for the ap- 
plication of generalised unitarity, [27] for a spacetime version thereof, [28, 29] for twistor 
techniques and [30] for the application of MHV diagrams. 

Computing form factors and correlators of non-protected local gauge-invariant oper- 
ators requires renormalisation, which in general implies the mixing of these operators. 
This procedure singles out certain subsectors, which are closed under renormalisation and 
which transform under subalgebras of the full PSU(2,2|4) symmetry [21]. The simplest 
testing ground for studying the full renormalisation problem of M = 4 SYM theory is 
given by the so-called SU(2) sector. The operators in this sector are built out of two 
complex scalar fields X and Y transforming in the fundamental representation of SU(2), 
eg. X = ġiı4 and Y = ¢gq4. In particular, the single-trace operators are of the form 
Obare = tr(X™Y*2 Xs y*4...), where kj € {0,1,2,...}. The renormalised operators of 
the interacting theory are obtained from these bare operators via the mixing matrix Z as 


Oren = ZObare ; Za 1472 +gtZ® + O(g). (1.2) 


The study of this mixing problem has been of great importance for capturing the novel 
integrable structures appearing in planar M = 4 SYM theory at higher loop orders [31]. 
At one-loop order, the crucial observation introducing integrability to planar M = 4 SYM 
theory was that the anomalous dilatation operator defined as 


co 
ôD = -n$ log Z = deg log Z = D (1.3) 
takes the form of the integrable Heisenberg spin-chain Hamiltonian within the SU(2) sector 
[32].? The central role of the dilatation operator and its interpretation as an (asymptotic) 
spin-chain Hamiltonian was further emphasized in [33], where the two-loop dilatation op- 
erator with SU(2) symmetry was computed from Feynman diagrams and its three-loop 
correction was derived under the assumption of integrability. A field-theoretic computa- 
tion of the latter was later performed in [34]. Making use of integrability, a recursive 
construction for the asymptotic dilatation operator in the SU(2) sector is available by 
now, which allows to compute its operatorial form to high orders in the ’t Hooft coupling 
constant [35, 36]. 
In this paper, we continue the program of [16, 17] and study form factors and the 
dilatation operator at two-loop order in the full SU(2) sector. We employ the unitarity 
‘Moreover, in [22] symmetry was used to show that all tree-level scattering amplitudes are related to 
certain contributions to the dilatation operator. The picture of [22] is equivalent to taking cuts of form 


factors. 
?In [32], the larger SO(6) sector was actually considered. 


method to obtain the complete two-loop form factors in this sector of planar VW = 4 SYM 
theory. Interestingly, the form factor results satisfy linear relations. It turns out that they 
can be explained by Ward identities of form factors following from R-symmetry. 

Form factors of non-protected operators contain both infrared (IR) divergences, due to 
soft and collinear virtual momenta, and ultraviolet (UV) divergences. The information of 
the latter allows us to determine the renormalisation matrix Z, and therefore, the dilatation 
operator.’ In dimensional regularisation, where the four-dimensional theory is continued 
to D = 4 — 2e dimensions,‘ all divergences are given by 1/e* terms. In order to obtain the 
dilatation operator, we need to disentangle the IR and UV divergences, which is possible 
since the IR divergences have a well-understood universal structure [37—40]. Concretely, 
we will subtract the IR divergences via the BDS ansatz [41, 42];° a similar procedure has 
already been used in [16, 17]. 

For amplitudes, it is well-known that the BDS ansatz does in general not give the 
full result but allows for a finite remainder function [18], which was first studied for the 
six-gluon case in [45-47]. For form factors of BPS operators, remainder functions have 
also been studied in [8, 14]. In particular, interesting properties associated to the so-called 
transcendentality were observed, such as the maximal transcendentality principle, which 
we will review below. In this paper, we will study the remainder functions of form factors 
of non-protected operators, where new features appear. 

Quantities in M = 4 SYM theory have shown interesting properties with respect 
to their transcendentality. Scattering amplitudes and form factors of BPS operators as 
well as their remainders have uniform transcendentality:® at ¢-loop order, they can be 
expressed as linear combinations of functions and numbers with transcendentality degree 
2L. Furthermore, remarkable relations have been found between the results of N = 4 SYM 
theory and QCD. It was first argued in [51] that, for anomalous dimensions of twist-two 
operators, the M = 4 SYM theory result is given by the leading transcendental part of 
the QCD result. This is usually referred to as the maximal transcendentality principle; 
see also [7, 52-54] for further discussions. While this heuristic relation was observed only 
for anomalous dimensions, in [8] it was found that the remainder function of certain BPS 
two-loop form factors matches exactly the leading transcendental part of related two-loop 
Higgs-to-gluons amplitudes in QCD [55]. This provides a first example where the maximal 
transcendentality principle is extended from pure numbers to functions which may have 
non-trivial kinematic dependence.’ 

In this paper, we demonstrate that form factors of non-protected operators show new 

3The anomalous dimensions can then be obtained as eigenvalues of the dilatation operator. 

“When continuing the spacetime dimension, also the fields have to be continued to D = 4 — 2e. This 
leads to some important subtleties which have been analysed in detail in [17]. These subtleties are, however, 
absent in the SU(2) sector. 

>See also the previous studies of amplitudes in QCD [43, 44]. 

° This is true at least in the cases of lower points or lower loops. There are known examples of amplitudes 
at sufficient high points in M = 4 SYM theory which are not given by transcendental functions but elliptic 
functions [48-50]. 


TAn interesting correspondence between the transcendental functions of M = 4 SYM theory and QCD 
was also found for energy-energy correlations [56]. 


universality properties regarding their transcendentality. Since the considered operators 
break supersymmetry, the remainder functions are expected not to have the property of 
uniform transcendentality. However, we find that all contributions of maximal transcen- 
dentality are identical to the corresponding results of BPS form factors. This provides 
further evidence for the universality of the leading transcendental part, which furthermore 
has a non-trivial kinematic dependence. 

This paper is organised as follows. In section 2, we present results for tree-level and 
one-loop form factors in the SU(2) sector. This also serves to introduce our conventions 
and notation. Moreover, we calculate the minimal two-loop form factors of such operators. 
In section 3, we extract the two-loop dilatation operator and two-loop remainder function 
from these results. Section 4 contains our conclusions and outlook. We provide simplified 
expressions for six-point amplitudes appearing in the unitarity calculation in appendix A. 


2 Minimal form factors in the SU(2) sector 


2.1 Tree-level form factors 


In this subsection, we summarise some general facts about form factors and give explicit 
tree-level expressions that are required in the unitarity calculations of the subsequent sub- 
sections. 

In analogy to amplitudes, we can strip off the gauge-group dependence of the form 
factors by introducing colour-ordered form factors Fo: 


Fo(l,...,m:q) = `> tr[T® ... T] Fo(a(1),...,0(n);q) + multi-trace terms, 
aESn/Zn 
(2.1) 


where T° with a = 1,..., N2 — 1 are the generators of the gauge group SU(N.) and the 
sum is over all non-cyclic permutations. The multi-trace terms in (2.1) can start to appear 
at one-loop order but are suppressed in the planar limit, and will not be considered in this 


paper. 
We describe the external on-shell states using Nair’s M = 4 on-shell superfield [57]: 
A,B Ay ByC 
7 non EABCD' N 
(p,n) = 9+(p) +0" palp) +y taB) +r 0° +n nnn g-i), (2-2) 


where 74 are Graßmann variables that encode the flavour and helicity of the component 
particles, and A = 1,...,4 is the SU(4) R-symmetry index. In this formalism, we can 
combine form factors with different external fields into one super form factor. As we will 
see later, this also makes it easier to study the supersymmetry properties of the form 
factors. 

In this paper, we focus on form factors in the SU(2) sector. The corresponding single- 
trace operators involve two complex scalar fields with a common SU(4) index, which are 
chosen explicitly as X = $14 and Y = @¢o4. The tree-level minimal super form factor for 
the operator O = tr(XXYX---) with L = n fields, for instance, is simply given by 


L 
0 x i p 
FP... Lq) = 5(q — YO Aike) (ninindndrgndnin} --- + cyclic permutations) . (2.3) 
i=1 


In general, the colour-ordered minimal tree-level super form factors of any operator can be 
obtained from the operator’s oscillator representation by replacing the oscillators by spinor 
helicity variables and multiplying the result by the momentum-conserving delta function 
[16]. 

We also need the next-to-minimal tree-level form factors in the two-loop unitarity 
computation below, which contain one more external field than the minimal ones. They 
may be computed easily by Feynman diagrams, or obtained from the BPS form factor 
component expressions, see e.g. [13]. For convenience, we provide some explicit rules that 
are useful in practice. There are four different cases that can occur. In the first case, a g4 
can be inserted between two neighbouring positions 7 and i+ 1. This leads to the following 
replacement in the colour-ordered minimal tree-level form factor: 

B (ii+2) 


A „B,C „D A © -D 
Oi Mi MAM OG Tiatia t2 ae (2.4) 


In the second case, a g_ can be inserted at the same position, leading to 


A,B,C „D A,B [i i+2] 1 2 3 4 C D 
TE N Nae i S eT Girip ipa i iiaiai sees (2.5) 
In the third case, a dcp at position i is split into two anti-fermions wc and Yp. This leads 
to 


1 C.D DO\E oF 
Gitl) (ni M41 — Mi NeA)Ni+2M 42 phe (2.6) 


A „B „C „D,E „F A „B 
o MiMi Mmi MiyiMi+i i S t Miia 
In the fourth case, the dcp is split into two fermions W@ and y?’ with ecpoæp = 1, 
leading to 
A „B „C„D,E „F A,B __— S e E F 
Hamam Mi Miyiz: > TAM gy Po TD p Mio itani oE 
(2.7) 
where ñi, A = He ABCDNËNE në and the minus sign is related to the order of the 7’s. The 
complete next-to-minimal form factor is obtained by summing over all four replacements 
and all insertion points. 


2.2 One-loop form factors 


In this subsection, we consider the one-loop minimal form factors in the SU(2) sector and 
show how to obtain the one-loop dilatation operator from them. This also allows us to 
introduce our notation and some important concepts that are required for the two-loop 
case. The results for the one-loop form factors, as well as the recipe to obtain the one-loop 
dilatation operator, were already given in [16]. Here, a useful new formulation, given in 
(2.11), is developed, which will be convenient to study the symmetry properties of form 
factors. 
Form factors in the loop expansion can be written in the following form: 


Fo = (1+ PIP 4 pTO +...) FY. (2.8) 


Pı 


P2 
Figure 1: The one-loop (pı + p2)? double cut. 


For operators that are eigenstates under renormalisation, such as BPS operators or the 
Konishi primary, Z is simply the ratio of the /-loop and tree-level form factor. However, 
for form factors of operators that renormalise non-diagonally, this is no longer the case, 
because the loop corrections to vanishing tree-level form factors can be non-vanishing. To 
overcome this problem, it is necessary to promote Z™ to an operator that acts on the 
tree-level form factor FS ) and creates a different tree-level form factor from it. 

In the planar limit, connected ¢-loop interactions can maximally involve Z + 1 neigh- 
bouring fields in the colour-ordered form factor at a time. Hence, Z can be written as 
an interaction density that is summed over all insertion points. At one-loop order, the 
maximal interaction range is two, and we can write 


L 
TS hore (2.9) 
i=1 


Here, L denotes the length of the operator O, r acts on the external fields 7 and i+ 1 


and cyclic identification i + L ~ 7 is understood. We depict w as 


Las @), (2.10) 


where we in general specify only the first field ¿i that is acted on when the range is explicitly 
specified by the number of occurring legs. 

In the SU(2) sector, the following six range-two interactions are allowed by R-charge 
conservation: XX > XX, XY > XY, XY > YX, YY > YY, YX > YX and 
YX — XY. It is sufficient to consider the first three, as the last three can be obtained from 
them by replacing X + Y, which is a symmetry of the theory. We denote the contribution 
to a given combination of external fields ZAZg > ZoZp by (zez, where Z1 = X, 
Z2 = Y and A,B,C,D = 1,2. In terms of these matrix elements, the operator Bs is 
explicitly given by 


(1) (1))\Z0Zp, C ð D ð 
I; I ~ —_ N —— , 2.11 
( zS zoni anA H ane, ( ) 


Table 1: Linear combinations of diagrams contributing to the minimal one-loop form 
factors in the SU(2) sector. 


The matrix elements (I (20% can be computed via unitarity. In the one-loop case, 
we only need to consider the double cut shown in figure 1. Let us briefly consider the 


(yx case. The cut integrand is given by 


f arsa, lo) d'm dm, FO EE ps,..- pr; QAO (np p), (2.12) 


where the tree-level form factor is given in (2.3) and the four-point amplitude is given by 
the standard MHV expression. The labelling of the external legs with X, Y in the tree-level 
amplitude and form factor means to take the corresponding 7 components; for example, 
A4(—l2,—l1, py , px) means to take the component of A4(—l2, —l1, p1, p2) containing the 
(n?nt)(ndnz) factor. Integrating out the m, variables, the cut integrand is given by® 


(aint bn FQ OT pa sPLADlydonger f APS). e 


The variables 7; and ņ2 indicate that the result is not necessarily proportional to the 
tree-level form factor of the original operator but to the one of the operator in which the 
corresponding X and Y fields are permuted. The occurring phase space integral is simply 


<x (2.15) 


At one-loop level, this cut is sufficient to determine the matrix element xx as the 


the cut of a scalar bubble integral: 


bubble integral. The other matrix elements can be obtained in a similar way. More details 
of such computations can be found e.g. in [16, 17]. 
The one-loop results are summarised in table 1. It is interesting to note that 


(LOVEE + DYE = PY. ee) 


v KA 
®Note that 


ə? a? 


(0); X Y g = (0) : 
Fo (pi > P2 P3» -o PL3 d) | Ann Ant = Oni Ont ono” © (pi,---,PL3q) nA=nA=o- 


This relation is a consequence of the SU(2) symmetry of the theory. Let us establish a 
formalism to deal with these symmetries in more detail since it demonstrates the general 
principle of how symmetries can be used to study form factors. 

The PSU(2, 2|4) symmetry of MN = 4 SYM theory leads to the following Ward identity 
of form factors: 


S o Foes) = Fzao(l,... nq), (2.17) 
i=1 


which holds for any generator J4 of PSU(2, 2|4); see e.g. [4] for a detailed derivation. Let 
us consider explicitly the generators’ 


o o o o 

~l 1 2 ~2 . 1 -2 ~3 1 2 

g = ~~ + E Ta = M; = tH iN —, 2 r a  _— nt —_ 218 
i i One i Onl Ji i One h Onl Ji i Onl M One ( ) 


of SU(2). Applying (2.17) to (2.8) for the minimal tree-level and one-loop form factor, we 
find 
krae (2.19) 


where J4 = Da 3A. Inserting (2.11) into (2.19) yields (2.16) as well as similar identities. 
The results of table 1 contain the one-mass triangle and bubble integral, for which 
explicit expressions can be found e.g. in [59]. The one-mass triangle integral is IR divergent 
and UV finite. The bubble integral, on the other hand, is IR finite but UV divergent. Hence, 
the IR and UV divergences can be separated immediately. 
The IR divergences of the above results match the universal form of one-loop IR 
divergences [1]: 


1 S 
Da R ~ a (—sii+1) € Lisi t+ O(e°) 
Ly g@ i (2.20) 
— cus se 
ne a 4e (—sii1) € ligt Ole); 


where yh, = 8 is the one-loop cusp anomalous dimension and gP = 0 is the one-loop 


collinear anomalous dimension. We have also introduced the identity operator 


Lii = 3 E (2.21) 
tt = i A i+ B : 3 
A,B=1 On; Ons 


The UV divergences require the renormalisation of the operators. The renormalised 
operators are defined in terms of the bare operators and the renormalisation constant Z 
as shown in (1.2). The renormalised form factor is nothing but the form factor of the 
renormalised operator.'? Since the form factor is linear in the operator, we can write in 
the case of the minimal form factor: 


Pome hii hesh (2.22) 


°In general, the generators of PSU(2,2|4) may obtain anomaly contributions, see e.g. [58]. These are, 
however, absent for SU(2). 

This statement relies on the finiteness of M = 4 SYM theory and on a formulation in which also 
wave-function renormalisation is absent. 


where, on the right hand side, Z acts as an operator on the tree-level form factor, similar 
to Z discussed before, cf. (2.8). 
At one-loop level, ZC) has to render the renormalised one-loop interaction 


TY = 7M 4 ZO) (2.23) 


UV finite. This means that z0 


occurring in J ee The UV divergence of the bubble integral is given by = Accordingly, 


has to cancel the UV divergence of the bubble integrals 


using the results in table 1, the one-loop renormalisation constant density is given by the 
matrix elements 


1 1 
(2% =0, ZDR PEs- (2.24) 


It can be written in the compact operatorial form 


1 
Ziti = = a(t =P iia (2.25) 


where 1 is the identity operator (2.21) and 


o 
Pitti = D ne (2.26) 


aa 4a B 
A,B=1 On ‘ On 


denotes the permutation operator. 

In analogy to the renormalisation constant, we can also write the dilatation operator 
as an operator acting on the minimal tree-level form factor. Applying (1.3) to (2.24), we 
find the one-loop dilatation operator density 


7 a 


These expressions can be combined into the well-known form [32] 
Dia = (1 Pig. (2.28) 
Let us now proceed to two-loop order. 


2.3 Two-loop form factors 


In the two-loop case, the range of connected interactions can be either two or three. Fur- 
thermore, two disconnected one-loop interactions can occur at two-loop level. In total, we 
can introduce the two-loop operator Z) similar to the one-loop case as 


L 
2) (2) 1) 1) 
T® = >, (ieat Ha +5 5 5 Ra Pal (2.29) 
i=1 F= =i+2 


"We use a modified minimal subtraction scheme with effective planar coupling constant g? = 


(4r e7)" Sae, 


10 


Sii+1Sil 


Table 2: Linear combinations of diagrams of range two contributing to the minimal two- 
loop form factors in the SU(2) sector. Terms between horizontal lines always occur in fixed 
combinations. 


where the last term accounts for the insertion of two one-loop interactions Ii at non- 
overlapping positions. The two-loop interactions 1), and D +2 are given by 
= (2) o O 
2 2)\Z02Zp.C D 
Ii = = `> (L; zS zpi gA HAB 
A,B,C,D=1 "i +A 
i 5 r 5 (2.30) 
(2) m _ (2)\Zp2ZnZr,D E F 
Lie = = > UI gaH BB 429,07: 

A,B,C,D,E,F=1 ni N41 N42 


For interaction range two, three distinct cases occur: XX > XX, XY — XY and 
XY — YX. For interaction range three, six distinct cases occur: XXX —> XXX, 
XXY —> XXY, XYX —> XYX, XXY > XYX, XYX > XXY and XXY => YXX. 
The remaining combinations can be obtained from these cases by exchanging X © Y 
and by using parity, i.e. reverting the order of the fields. We collect our results for the 


1))xx 


corresponding matrix elements in table 2 and table 3. The matrix elements (J;°’)<{% and 


(P )xxx occur in the BPS case and were computed in [14] using the unitarity method. 
The other matrix elements can be calculated in a similar but slightly more involved com- 
putation. Let us give a brief account of this computation. 

The cuts that have to be considered are depicted in figures 2a, 2b, 2c and 2d. The tree- 
level next-to-minimal and one-loop minimal form factors, which occur as building blocks, 
are given in subsections 2.1 and 2.2, respectively. The required tree-level and one-loop 


amplitudes are standard. A particularly interesting cut is the triple cut shown in figure 2d, 


11 


XXX XXY AYA XYX XXY 
XXX XXY AYA XXY XYX 


+1 
+1 
-1 
+1 
+1 
-1 
+1 
+1 
+1 
-1 
1 


SilSi+li+2 | +1 | +1 
Sii+1Si+21 | +1 | +1 
Sii+li+2 -1 -1 
+1 
+1 
-1 


DERE 


Table 3: Linear combinations of integrals of range three contributing to the minimal two- 


Si+1i+2 


loop form factors in the SU(2) sector. The integrals are grouped such that those between 
two horizontal lines are always occurring in the same combination. The second and fourth 
group as well as the third and fifth group are related by parity. 


which involves the tree-level next-to-MHV six-point scalar amplitudes. As given explicitly 
in appendix A, these scalar amplitudes take a simple form in terms of Mandelstam variables. 


[xxx 


Let us consider for example the (J;°’)\+ case. The cut integrand is given by the product 


12 


(c) The two-loop (pı + p2)? triple cut. (d) The two-loop (pi + p2 + p3)? triple cut. 


Figure 2: Unitary cuts of the minimal two-loop form factor. 


of the minimal tree-level form factor and the six-point tree-level amplitude: 


3 
[cutest tot) [ft FS 5 Fa, ..., pti DAg(—ls, —l2, —h, pi, p% ,pž). 


i=1 
(2.31) 
After integrating out the 7, variables, we obtain 
0 
FS) (pit PY PS +--+ sPLD)| yA py (MME) (TMB) (NBM) 
1 (2.32) 


x farsa, lls) 


where we have used (A.2) for (A)***.!2 The phase space integral corresponds to the 


ANY to, (2.33) 


l3 
Similarly, each term in the other amplitudes in (A.2) is mapped to one graph in table 3 


triple-cut loop integral 


via the triple cut. 
Looking at table 2 and table 3 and using the parity transformation, we observe some 
linear identities, e.g. 


(1 yxx UP =, (2.34) 
and 2) 2) 2) 2) 
C30 ed ar 
oan T a cee. OH a neem (2.35) 


12Note that there is a minus sign from the m, integration which cancels the sign in (A()¥X*. 
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Similar to the one-loop case, all these identities are a consequence of SU(2) invariance and 
follow from the Ward identity (2.17), which at two-loop order yields 


ot 0: (2.36) 


Given the full integrand of the two-loop form factor, we can perform a similar analysis 
as in the one-loop case. However, we will see that this requires a more involved subtraction 
of the IR divergences. This will be the topic of the next section. 


3 Two-loop dilatation operator and remainder function 


In the one-loop case, the UV divergences stem from the bubble integrals alone. Therefore, 
the one-loop renormalisation constant can be read off directly from the coefficient of these 
integrals. This is no longer true for two-loop form factors, since the two-loop integrals 
in general contain a mixing of IR and UV divergences. However, IR divergences have a 
well-understood universal structure [37-40]. This allows us to subtract the IR divergences 
systematically using the BDS ansatz [41, 42]. 

Similar to the one-loop case (2.23), the two-loop renormalised form factor is given by 


T? —72 +720 4 ZO, (3.1) 
where 7 
1) (1) 
Ze?) = 5 C ii+li+2 +5 9 5 20,20). (3.2) 
i=1 32 =i+2 


Applying the BDS ansatz [41, 42] to the renormalised form factors, we obtain a finite 


two-loop remainder function: 
2 
RÐ = IP (e) — 5 (T90) - SPE) + Of), (3.3) 


where 
fË (e) = 28 — Lze = Xe? . (3.4) 
At two-loop order, connected interactions involve at most three fields of the composite 
operator, which have to be adjacent at the planar level. Hence, both the remainder function 
and the dilatation operator can be written in terms of densities that act only on triples of 
neighbouring sites at a time and are summed over all L insertion points. For each triple of 
neighbouring points, we define the variables 


Sii+l Si+li+2 Si+2i 
y= ys E, m= H, (3.5) 
Sii+li+2 Sii+li+2 Sii+li+2 
where 
Sii+li+2 = Siitl + Silit? + Si+2i (3.6) 


and cyclic identification į ~ i + L is understood. These variables satisfy u; + vi + w; = 1. 
The remainder RC) can be written in terms of its density as 


L 
or YR ie (3.7) 
i=1 
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An important subtlety arises due to the fact that the composite operators are not 
necessarily eigenstates under renormalisation. This requires a careful treatment of the 
product of one-loop form factors in (3.3). As already mentioned, the renormalisation 
constant is a matrix (i.e. an operator) and so are the interactions. Hence, the one-loop 
product in (3.3) should be understood as a product of operators. This can be explicitly 
depicted by the following equation in terms of graphs: 


L L+i-2 
wo -2 H POE BS OE) - (3.8) 


Note that the states corresponding to the internal lines are summed over as required for 


a product of operators. The prefactors of 2 stem from distributing products of densities 
with effective range two equally between the first two and the last two sites. 
The remainder density, which itself is an operator, can be similarly expressed by the 


4 mn TORON 
(3.9) 


—2e 
where we have depicted the two-loop renormalisation constant density Zz), j42 analogously 


following graph equation: 


2 
fies a 


to 1, 49. Requiring that this remainder density is finite allows us to fix the two-loop 
renormalisation constant density. 

The integrals occurring in the two-loop result can be reduced to master integrals via 
IBP reduction, e.g. as implemented in the Mathematica package LiteRed [60]. The result- 
ing master integrals can be found in [61]. 


3.1  Renormalisation constant and dilatation operator 


From the requirement that the two-loop renormalisation constant densities have to cancel 
all divergences in (3.9), we find 


2 2 2 2 2 1 1 
(2.9% = 0, (2k =t+a-s CPt a 
2 1 1 2 1 1 (2) 1 1 
Co oau E o a a S -Z- 
(3.10) 
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Alternatively, this can be written in the operatorial form! 


1/1 1 
ee =o (= = :) (PiitiPipiito + PipiitoPiigi — 3Piipi — 3Pipiit2 +4). (3.11) 
Using (1.3), we have for the two-loop dilatation operator 


D?) = de (2® = 5(2)) (3.12) 


where (Z (1))2 should be understood as a product of operators. For example, 


2 2 1, 1 1, 1 
OPi = de (a a l aa o 
(3.13) 
1 1 1 1 1 1 
-FEMRE - EDERE) =4. 
In total, we have 
2 2 2 
@PIE=0,  @PYEE=—-8, PVE =-2, gw 
2 2 2 i 
oers OPR =, OPS Ee 
which agrees exactly with the known result [33] 
D a = —2(Pii p1Pi+ii+2 + PigiioPiigi — 8Piin1 — 3Pi+ii+2 + 4) ; (3.15) 


3.2 Finite remainders 


Next, we calculate the finite remainder densities. The remainder densities fulfil analogous 
relations to (2.35): 


(RE P VERE + (RP VY = (RPV, 
RË XYX RË YXX RË XXY _ RË XXX (3.16) 
( 4 are | 4 ee i \xyx = ( i Jex: > 


aR So, (3.17) 


which is a consequence of (2.19) and (2.36). Combining (3.16) with the symmetry under 
the exchange of X + Y and the reversion of the order of the fields, we can express all 
remainder densities in terms of (RO))XXX, (RO) x¥x and (RP XXX, Hence, it is enough 
to consider these three cases. 


13Note that the coefficient of the simple pole coincides with the one of the double pole up to a sign. 


This is a consequence of the fact that at two loops only one Feynman integral with overall UV divergence 
occurs in the SU(2) sector when using a manifestly IR finite formulation, and it does not hold for general 
operators; cf. [34]. 
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The remainder density (RO )Xxx was already studied in [14] and is of homogeneous 


transcendentality four: 


2 2 
(PG as (3.18) 
which is given explicitly as 
2 
Gre, 
gol 1— ui . gol . 
= —Li4(1 — Ui) — Lis (u;i) + Lig ( ) — log ( “ ) [ti (- ) — Lig (1 — w) 


i : i ‘ i— l1 1 1 
= log (ui) [tis (=) + Lis (-=) + Lis (: D ) — 3 log? (vi) = 3 log? (1 = w) 


uj; —1 Vi 1 — ùi 1 1-4; 
=i 4 Li i Li ; l i l : L] 2 i 
i ( m ) ip (_)+ in (us) [log (= ) og (vi) + 5 og (= )| 


1 1 1 A 
+5 low! (us) — G log? (u) 10g? (wi) — 5 10g? (1 — us) og (u) 10g (=) 


7 


1 1 
=z log (1 — us) log” (uj) log (vi) — 7 log? (uj) log (w;) 


—@| log (ui) log (- = *) + 5 log? (- —*) = 5 log? (ui) | 


+¢3 log(u;) + “ + G({1 — uj, 1 — ui, 1,0}, vi) + (ui © vi). (3.19) 


Here, the Goncharov polylogarithm in the last line is the only piece that cannot be written 
in terms of classical polylogarithms. This relatively compact expression was obtained using 
the symbol techniques [62, 63]. The corresponding symbol is given by [14]'4 


uj-1l_ uv Ui 2 l-u; w 
S ((RP Ex |) = -u De i Vi i il i Wi 
(Ri )xxXx|, uj @ (1 — ui) 8 : ae” a es ui 8 Uj ® are 
Ui Uy Ui Ui 
-ui Q vi ®@ — 9 — — u 8 vi D — 8 — + (u; ou). (3.20) 
i Wi i Wi 


The remainder density (RO? )x¥x is of mixed transcendentality with degree ranging 


from three to zero. Its contribution of degree three has the symbol 


vi Ug l— vi _ vi l— ui _ vi 
S (RP aal) = ; t t o tg = % *@—. (3.21 
(Ri )xxY | Taa TE vi 8 Ui E ue Ui a ( ) 


The full transcendentality-three part can be given as 
. Uj f vi 1 we 
(APs = [is (-4) ~ Jog (ui) Lis (=) + log (1 — wi) log (u) log G i ) 


1 iUi 1 1 
— z (-= ) E log (u;) log (vi) log (w;) — — log? (w;) + (u; © v| 


l 


1 l— vi 1 l 
-Liz (1 — v;) + Lis (u;) — = log? (v;) log (=) + <1? log (=) 
2 Uj w 
1 
3a log (—si i+1 i+2) ‘ (3.22) 


14The symbols can be conveniently calculated using the Mathematica code [64]. 
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Together with the terms of lower transcendentality, we have 


7 (3 


RY xxx — (R@ xxx) 4 Lig (1 — uy) + ilii (3.23) 
XxX 3 


2 


1 Uj T 
+ log (ui) log (vi) = 5 log (—Si+1i+2) log (=) + 2log ( Siiri) } 3 T. 


The final remainder density (RP xxx is of mixed transcendentality with degree rang- 


ing from two to zero. It reads 


1 Ui 4 1 
(RPE = Flog (-sisasea) tog (SE) — Lio (1 — ui) — tog (ui) 10g (vi) + $ 10g? (0) 


T 
+ log (—5i+1i+2) — 2log (—sii+1) + 3° (3.24) 


Let us emphasise that, if non-vanishing, the transcendentality-four contribution is the 
same for all remainder function densities. Furthermore, there is only one transcendentality- 
three function and two functions of transcendentality smaller or equal to two that contribute 
to the results in the SU(2) sector. Notably, the highest degree of transcendentality t = 4— s 
is directly related to the shuffling number s of the respective remainder density, i.e. to the 
number indicating by how many legs the field flavours are shuffled. For instance, (RO? )x¥x 
has shuffling number s = 1 and maximal transcendentality degree t = 3. 

Interestingly, the rational pieces of the remainder function are connected to the dilata- 


tion operator as 


4 


(2) — (2) 
Diit = Zy Piili : 


: (3.25) 


4 Conclusion and outlook 


In this paper, we have calculated the two-loop minimal form factor for all operators in the 
SU(2) sector of planar M = 4 SYM theory via the on-shell method of unitarity. More- 
over, we have extracted the corresponding two-loop remainder function and the two-loop 
dilatation operator from it. The results of this paper provide a solid stepping stone towards 
calculating the complete two-loop dilatation operator of N = 4 SYM theory. The employed 
method, however, is independent of the high symmetry of planar M = 4 SYM theory, in 
particular of its integrability, and it is thus also applicable to less symmetric theories. 

The SU(2) sector is the simplest closed sector of the theory whose operators do not 
renormalise diagonally. It is hence well suited to study the occurrence of operator mix- 
ing and the dilatation operator. Due to the on-shell nature of the external fields, the 
divergences of the form factors are a combination of UV and IR divergences. We have 
disentangled the UV divergences from the IR divergences using the BDS ansatz and the 
universality of the latter. Because of the operator mixing, we needed to promote the inter- 
actions to operators and the iterative structure of the BDS ansatz to an operatorial form 
as well. From the UV divergences, we have determined the renormalisation constants and 
the dilatation operator. 

In contrast to the BPS case, the two-loop remainders of non-protected operators in the 
SU(2) sector do not exhibit maximal uniform transcendentality. However, their maximally 
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transcendental part coincides with the remainder of the BPS vacuum computed in [14]. 
These results and further evidence in other sectors lead us to conjecture that the two- 
loop remainder of every minimal form factor has the same degree-four part as the BPS 
one. Moreover, the two-loop remainder of the three-point form factor of every length-two 
operator should agree with the corresponding BPS remainder found in [8]. It would be 
interesting to check our conjecture about this universality for a wider class of operators, 
or even to prove it.!° 

Another observation is that the maximal transcendentality of the various remainders 
is related to their shuffling number, i.e. to the number indicating by how many places the 
field flavours are shuffled. It would be interesting to explore this pattern in larger sectors 
and at higher loops, where more complicated interactions contribute. 

Soft or collinear limits of scattering amplitudes or form factors are typically given 
by lower-point amplitudes or form factors multiplied by a universal function. Thus, they 
provide important constraints and sometimes even allow to bootstrap the full function 
under consideration. As already observed in [14] for the BPS case, soft and collinear limits 
of minimal form factors also do not vanish for the cases considered here. This may be 
surprising since a priori there is no physical interpretation for this limit because these form 
factors correspond to the minimal physical configuration. Actually, similar questions also 
appear at the amplitude level. For example, taking certain soft or collinear limits of the six- 
scalar amplitudes given in appendix A does not generate any physical amplitude, though 
the limit is non-zero. Via unitarity cuts, this consideration on the level of amplitudes 
affects the limits of minimal form factors. It would be interesting to understand this point 
better, and to see if one can obtain the soft or collinear limit without computing the full 
quantity. 

We have seen that the two-loop form factors, as well as the remainder functions and 
dilatation operator derived from them, obey Ward identities induced by the underlying 
R-symmetry. Going beyond the SU(2) sector, the realisation of similar Ward identities fol- 
lowing from other symmetries should be more involved due to corrections to the generators 
which are absent for SU(2), cf. examples of such corrections in the case of spin chains [66] 
or scattering amplitudes [58, 67]. Capturing these extensions in the case of form factors is 
currently under investigation. 
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A Six-point scalar amplitudes 


In this appendix, we provide all six-point amplitudes that are required in the unitarity 
computation of the two-loop form factors in the SU(2) sector. We use the short notation 


(AO zzz = AO (171,222,323, 474,525,625) a> (A-1) 


where, on the right-hand side, all momenta are taken to be outgoing and we have set all 
ņi variables to 1. 
The required amplitudes can be explicitly given in terms of Mandelstam variables as 


(AOR =- —, 
aog = — —, 
(A@)z2% = 523556 4 512545 8123 l 
5168345234 $16 5345345 516634 
(aO = a OL L, 
8168345 8168234 $16 $8345 
(AM = 4 8 
5345234 8345345 $34 $234 
(AM xa¥ Fe _*50 eee TE (A.2) 
6165345 S165234 $16 $234 
TO =—B_ 4 Sg 
5345234  §345345 534 5345 
(A))xx¥ = 523856 512945 4 5123 $12 S56 aa 
8169345234 5165345345 $16834 5165345 $165234 S16 
(A@yexx $s Ssa $128 __—823_ 8a E 
§165345234 5165345345 $16534 8349234 $345345 534 
(AO) x¥x = — 523556 912545 4 $123 
§165345234 5165345345  $16534 
a a E E SSE S 
5168345  $165234 $16 $345 5345234  $345345 S34 $234 


where all poles are physical.!° Through the triple cut shown in figure 2d, each term in 
(A.2) is exactly mapped to one graph in table 3. 


‘©The BCFW recursion relation or MHV rule methods directly give results in a much more complicated 
form, usually involving spurious poles. We have checked that these different methods give results equivalent 
to (A.2). 
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It is easy to find various relations among these amplitudes, such as 


(O0)\y¥xx |. (0)\xyx (O)\xxy _ _/ 4(0)\Xxx 
(AN) (AN) (A) (AN) 


XXY F 


XXY XKY ~ XXX) 
(AO a A ua = (A.3) 
(A ee HAM = (Aa ee 


These are the counterparts of (2.35) and nothing but supersymmetric Ward identities 
(SWI) for amplitudes [68].17 
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